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Abstract 



It is shown that the law of an SDE driven by fractional Brownian motion with Hurst 
parameter greater than 1 jl has a smooth density with respect to Lebesgue measure, provided 
that the driving vector fields satisfy Hormander's condition. The main new ingredient of the 
proof is an extension of Norris' lemma to this situation. 
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1 Introduction 

In the celebrated paper IHor67l . L. Hormander gave in 1967 a sufficient (and necessary in the 
analytic case) condition for the hypoellipticity of second order differential operators. The original 
proof of Hormander was rather complicated and has been since considerably simplified in using 
the theory of pseudo-differential operators. 

In 1976, P. Malliavin, uses the deep connection between the theory of second order differential 
operators and Ito's theory of stochastic differential equations to point out the probabilistic coun- 
terpart of Hormander's theorem. The problem of the hypoellipticity is essentially equivalent to the 



1 



problem of the existence of a smooth density with respect to the Lebesgue measure for the solution 
of a stochastic differential equation. The idea of Malliavin's proof of Hormander's theorem is to 
show that the Ito's map associated with a stochastic differential equation is differentiable in a weak 
sense and then to show that, under Hormander's conditions, this derivative is non-degenerate. The 
stochastic calculus of variations that has been developed in IMal78l precisely in order to obtain 
a probabilistic proof of Hormander's theorem is now known as Malliavin calculus and has since 
then found numerous applications (see INua95IO . 

In the last few years, there have been numerous attempts to define a notion of solution for 
differential equations driven by a fractional Brownian motion. When the Hurst parameter of the 
fractional Brownian motion is greater than 1/2, existence and uniqueness of the solution are ob- 
tained by Zahle in HZahOll or Nualart-Ra§canu in IINR02I . Let us note that, as a consequence of 
the work of Coutin and Qian |CQ02| , a notion of solution can actually be well-defined for H > |. 
The problem of the existence and smoothness of the density with respect to Lebesgue measure for 
solutions of stochastic differential equations that are driven by a fractional Brownian motion with 
Hurst parameter greater than 1/2 is solved in some special cases. In [NT06] the existence and 
smoothness of the density has been shown in the one-dimensional case by using Doss-Siissman 
methods. In INS05I . the authors prove the existence of a density under ellipticity assumptions. Fi- 
nally, in IHN06I . always under an ellipticity assumption, the smoothness of the density is proved. 

In the present work we prove a version of Hormander's theorem for solutions of differential 
equations driven by a fractional Brownian motion. More precisely, we prove that as soon as the 
usual Hormander conditions are satisfied, there exists a smooth density for the law of the solutions. 
One of the major difficulties is to prove an equivalent for fractional Brownian motion of Norris' 
lemma INor86l . This lemma, interesting in itself, quantifies in which way the integrand of a 
stochastic integral with respect to fractional Brownian motion has to be small if the stochastic 
integral is small. 

This article is organised as follows. In a first section we recall the basics of Malliavin calculus 
in the context of fractional Brownian motion, which is the main tool used in this article. We 
then prove a version of Norris' lemma for stochastic integrals with respect to fractional Brownian 
motion. This lemma is then used to prove the main result of this work, that is the existence and 
the smoothness of the density under Hormander's conditions. Finally, in a last section, we apply 
this result to the analysis of the behaviour in small times of this density on the diagonal. 
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2 Malliavin calculus with respect to fractional Brownian motion 

Let us first recall some basic facts about Malliavin calculus with respect to the fractional Brownian 
motion (for further details, we refer for instance to IHN06I or [NS05|). 
We consider the Wiener space of continuous paths: 
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where: 

1. C([0, 1], R d ) is the space of continuous functions [0, 1] ->■ R d ; 

2. (A) t>0 is the coordinate process defined by /3 t (f) = f(t), f G C([0, 1], R d ); 

3. P is the Wiener measure; 

4. (Bt)o<t<i is the (P-completed) natural filtration of (A)o<t<i- 

A ci-dimensional fractional Brownian motion with Hurst parameter H e (0, 1) is a Gaussian 
process 

B t = (Bl...,Bf), t>0, 
where B l , . . . , B d are d independent centred Gaussian processes with covariance function 

R{t ,s)= l -(s 2H + t 2H -\t-s\ 2H ). 

It can be shown that such a process admits a continuous version whose paths are Holder p contin- 
uous, p < H. The fractional Brownian motion with Hurst parameter H > ~ can be constructed 
on the Wiener space by a Volterra type representation. Namely, the process 

B t = [ K(t,s)dp s ,t>0 (2.1) 
J o 

is a fractional Brownian motion with Hurst parameter H, where 



(u- s) H ~^u H ~^du , t 



> s. 



and ch is a suitable constant. Let £ be the space of M. d - valued step functions on [0, 1]. We denote 
by TL the closure of £ for the scalar product: 

d 

((^•[o.ti]) • • • i ^-[o.td])) (ip.siji • • • A[Q, Sd ]))n = y~]R{ti,Sj). 

i=i 

For tp,ip eH, we have 

V)« = ff(2i? - 1) / / | s - t \ 2H ~ 2 (<p(s)M*))jL*d*&- 
Jo Jo 

It can be shown that L 1 /^([0, 1], M d ) C TL but that TL also contains distributions (see for example 
IpToTH ). 

For our purposes, the following representation of the H-scalai - product is useful. Let T a ip 
denote the fractional integral of order a of tp, defined by 

l a if(t) = -i- f (t - s) a - V(s) ds . (2.2) 
"(a) Jo 
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If one extends ip and ip to K+ by setting ip(t) = tp(t) = for t > 1, a straightforward application 
of Fubini's theorem shows that one has the identity 

{^) n = (l H - l /\,l H ~ 1/2 i>) L > • (2-3) 
A Si -measurable real valued random variable F is said to be cylindrical if it can be written as 



F = f([ (hl,dB s ),..., [ (h^,dB s 



where h l G Ti and / : R n — ► M is a C°° bounded function. The set of cylindrical random variables 
is denoted S. The Malliavin derivative of F G 5 is the R d valued stochastic process (D^F)o<i<i 
given by 

n t F = hi ( t )§^.{J Q < h ^ dB '^ ■■■J o (K,dB s 

More generally, we can introduce iterated derivatives. If F G 5, we set 

T>l_ tk F = T> tl ...T> tk F 

For any p > 1, the operator D fc is closable from <S into L p (C([0, l],R d ), We denote by 

J} k,p (TL) the closure of the class of cylindrical random variables with respect to the norm 

\\n k , P = Un+^E(||D^n^, 

and 

b°° (h) = pi pi B fc ' p (w). 
p>i fe>i 

We also introduce the localised spaces D^(7i) by saying that a random variable F belongs to 
Dj (7i) if there exists a sequence of sets Q n C B\ and random variables F n G B fe,p (W) such that 
O n t C([0, 1], ]R d ) almost surely and such that F = F n on Q n . 

We then have the following key result which stems from Theorem 2.1.2 and Corollary 2.1.2. in 
INua95l : 




Theorem 2.1 Let F = (Fi, . . . , F n ) be a B\-measurable random vector such that: 

1. For every i = 1, . . . , n, Fi G 

2. 77ze matrix T = ((DF l , Di 7 ^)-^) 1<4 • <n w invertible almost surely. 

Then the law of F has a density with respect to the Lebesgue measure on W 1 . If moreover F G 
B°° (H) and, for every p > 1, 

^/n's density is smooth. 
Remark 2.2 The matrix T is called the Malliavin matrix of the random vector F. 
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3 Norris' lemma for integrals with respect to fractional Brownian 
motion 



A main ingredient in many probabilistic proofs of Hormander's theorem (see for example INua95IO 
is Norris' lemma [Nor86|. Loosely speaking, it is a more quantitative version of the uniqueness 
property of the semimartingale decomposition, stating that if a semimartingale is small, then both 
its bounded variation part and its martingale part must be small. In other words, the martingale 
part and the bounded variation part cannot compensate each other. This section is devoted to the 
proof of Proposition l3.4l which is a version of Norris' lemma formulated in a framework suitable 
for the purposes of this article. 

We start by stating a minor variant of a well-known concentration result for Gaussian measures. 



Lemma 3.1 Let fibe a Gaussian measure on a separable Hilbert space TL with covariance oper- 
ator r and write T = \ZtrT and A = \/||r||. Then, there exists a constant C independent of V 
such that one has the bound 



li{\\\x\\-T\>h) <Cexp{-^) 



for every h > 0. 



Proof. We can assume TL to be finite-dimensional; the general case follows from a simple ap- 
proximation argument since none of the constants depends on the dimension of TL. Define T = 
f \\x\\ p,(dx). It follows from the isoperimetric inequality for Gaussian measures ICIS761 ITal95l 
that 

li{\\\x\\-T\ >h) <Cexp(-^2j , (3.4) 
for every h G K. Furthermore, it follows from |Bog98 Thm 1.7.1] that 

r^o r^o VT 2 >9 , , m ~ 7T 2 A 2 7T 2 A 

T 2 - T 2 < —A 2 , and thus T -T < < . 

- 4 - 4T - 4 

The claim follows at once. □ 

We now define the class of Gaussian processes that are of interest to us. For H € (|, 1], we say 
that a Gaussian process B is of type H if it is centred and the function / defined by 

f(s,t) = E(B(t)-B(s)) 2 

is C 1 and satisfies 

ci\t-s\ 2H < f(s,t) <c 2 \t-s\ 2H , \d s dJ(s,t)\<c 3 \t-s\ 2H - 2 , (3.5) 



for every pair of times s,t £ (0, 1) with s ^ t. The following is a direct consequence of (13.5 
combined with 13. II 
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Lemma 3.2 Let B l be i.i.d. Gaussian processes of type H > ^ and let 5,N > be such that 
5N < 1. Define the W N -valued random variables X 1 and the number T by 

N 

Xi = B\n5) - B'dn - 1)5) , T 2 = £ f(nd, (n - 1)5) . 

71=1 

77je7i, ?/iere e»'s? constants C\, C2 such that one has the bound 

( Nh 2 \ 

P(||X»| -T\ >h)< C.exp (-C 2 p^J, (3.6) 

Nh 2 
(6Nf 



\(X\X 3 )\ > h 2 ) < C ieX p ( -C 2 — m ), (3.7) 



for every h > and every pair (i, j). 

Proof. Denote by T the covariance of X 1 (this is independent of i). Then, one has 



r m ,n = - / / d s d t f(s,t)dsdt, 

^ J mS J n5 



so that |r mn | < C5 2H (1 + \m - n\) 2H - 2 . This implies that 

N N 

\\r\\ls= E \T mn \ 2 <CN5 m Y J \k\ AH ~* <C5 m N m ~ 2 . (3.8) 

m,n=l k=l 

Since this is a bound on ||r|| 2 , the first inequality follows from Lemma l3~T1 Fix now an arbitrary 
pair of indices i ^ j. Note that conditional on the value of X 3 , the random variable (X l ,X 3 ) 
is normal with variance (X 3 , FX 3 }. This motivates the introduction of the random vector X 3 = 
Y X I 2 X 3 which is Gaussian with covariance T 2 . Note also that EHX 7 '!) 2 = tr(T 2 ) = ||r||^ s . We 
thus have, for any v > 0, the bound 

P(\(X\X j )\ >h 2 ) =E(F(\(X\X j )\ >h 2 \X 3 )) <E(Cexp(-/i 4 /||^ i || 2 )) 

^ Cex P (- Jnv^nZ , ^2 ) + P (HI^II " H r llHs| > v) 

v 2 



h 4 


4(l|r| 


Ihs + v) 2 




h 4 


8(||r| 


Ihs + - 2 ) 




h 4 


8(l|r| 


Ihs + - 2 ) 



< Cexp ( - „ ^ + Cexp 



4||r 2 

< C exp ( - =- ) + C exp 



4||r||| s 



Note now that the second inequality is non-trivial only for h 2 > ||r ||hs» so that we assume that 
we are in this situation from now on. Choosing v 2 = /i 2 ||r||HS> we g et 

P{\(X\X 3 )\ > h 2 ) < Cexp (-/i 2 /(16||r|| H s)) (3.9) 
which, together with (I3.8I) . implies the required bound. □ 
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Remark 3.3 The bounds in Lemma 13.21 can be interpreted as saying that the 'coarse-grained 
quadratic variation' of B on a scale 5 and over a time interval t behaves like 5 2H ~ 1 t to within 
an error of order 5 H t H . Note that the relative magnitude of the error to the average value always 
tends to as 5 — > 0, but that this ratio becomes 'worse' as H — > 1. 

We now have the main tools in place to prove the following version of Norris' lemma. Note 
that here and in the sequel, we denote by || • || Q the a-Holder norm of a function. 

Proposition 3.4 Let H G 1) and let a and b be processes taking values in R and R m respec- 
tively such that'E(\\a\\g + ^ ||& l ||#) P < oo for every p > 1 and every H S (^,H). Let 

yt= [ a(s)ds+ [ (b(s),dB(s)) , 
Jo Jo 

where the Bi are m i.i.d. Gaussian process of type H. Then there exists q > such that, for every 
p > 0, the estimate 

F(\\y\\L°° < e and \\a\\ L oo + \\b\\ L oo > e q ) < C p e p (3.10) 
holds. The constant C p depends on a, b, and p but not on e. 

Remark 3.5 Note that we do not require the Bi to be independent of the processes a and b. We 
also do not require any adaptedness at this stage. The reason why we will require adaptedness 
later on is that equation (14.211) for the Malliavin derivative of the solution does not hold otherwise. 

Remark 3.6 The bound (I3.10t actually implies the bound 

F(lblU- < e) < C p e p + min{P(||a|| L oo < £ «),F(\\b\\ L oo < e")} , (3.11) 
which will be used repeatedly in the sequel. 

Proof of Proposition MAX The proof consists of two parts. In the first part, we show that one has a 
bound of the type 

P(||y|| L oc < e and ||6|| Loo > e q ) < C p e p . (3.12) 
In the second part, we use this information to show that one has also 

P(||y|| L oo < e and ||a|| L oo > e q ) < C p e p . (3.13) 

Combining both bounds then yields (13.10b . 

The key idea to the proof of J3.12I) . which should be translated to 'if y is small then b must also 
be small' is the following. Choose two small length scales i5 « A C 1. Since we assume some 
regularity on b, it is easy to control the error made by assuming that b is constant on intervals of 
length A. One then considers the square root of the coarse-grained quadratic variation of y on 
a scale 5 over an interval of size A around t. This is of course bounded by ||y||^oo(5 _1 / 2 A 1 / 2 . 
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By Remark 1331 the contribution of the term including b to this expression is approximately equal 
to 5 H ~ 1 ^ 2 A l / 2 \b{t)\. The contribution of the term including a on the other hand is bounded by 
Summing over all intervals of size A yields a bound of the type 



< 6 



-Hi 



v\\l°° + 5 1 H \\a\\L 



from which it is then straightforward to deduce ( 13.121 ) by making use of the a priori bounds on the 
Holder norms of a and b. 

This argument is of course extremely sloppy, since we have not justified in any way some of 
the approximations made and we have not addressed the fact that b takes values in W 71 . Fix some 
small value of e > and fix two small numbers 5 and A such that 1/5 and 1/A are integers, 1/A 
divides 1/5, and such that (S << A < 1. We also fix if £ to be determined later. We 

define b(t) as the stepfunction with steps of length A approximating b(t), b(t) = 6(A[i/A]) and 
we write /3(t) = b(t) — b(t). The stochastic integral of j3 is bounded as follows. 

Lemma 3.7 There exists a constant C such that, for every t > and every s G [0, A], one has 

+S ((3(r),dB(r)} < C||6||^||5||^aM . 

Proof. Since s < A, we can assume without loss of generality that b is constant on the interval 
[t, t + s]. It follows from IYou36l that, for every a > 1, there exists a constant C such that 

" S (P(r),dB(r)} - (p(t), B(t + s) - B(t)) < 



h\\ b \\h s 



2H 



Furthermore, one has \ {f3(t),B(t + s) - B(t))\ < \\b\\ fl \\B\\ fl A H s H , so that the result follows 
at once. □ 

Denote by r the (integer) ratio A/ 5, set t n = 5n and define for N = 1, . . . , A -1 and for 
i,j = 1, . . . , m the random variable 

Nr-l 

X N= (A(Wl) - Bi(tn)) (Bj(t n+1 ) - Bj(t n )) . 

n=(N-l)r 

Note that T 2 := EXjj = rf(5) w A^ 2 ^" 1 . With these notations and using Lemmal3~71 we have, 
for n G \{N — l)r, Nr — 1], the relation 

\(b(NA),B(t n+1 )-B(t n ))\ < \y tn+1 -ytJ + \\a\\ L ~S+~ 



(P(s),dB(s)) 

< 2\\y\\ L °° + \\a\\ L °°5 + C\\b\\ 6 \\B\\jjA H 5 H . 

Taking squares on both sides, summing from n = (N — l)r to n = Nr — 1 and taking square 
roots, we get 



l^2b i {NA)b j (NA)X% < A 1 / 2 r 1 / 2 (2||y|| LO o + \\a\\ L o a 5 + C\\b\\ 6 \\B\\ 6 A 6 d A ) . (3.14) 
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At this point, it is convenient to introduce quantities Yfo = (X^) 1 / 2 and Y$ = {X^l 1 / 2 . With 
this notation, we get 



Y,\k(NA)\Yh < CY,Y^MNA)b,(NA)\ + CA 1 / 2 6- 1 ' 2 \\y\\ L «> 

i=l i^j 

+ CAVWWaUo* + C7A^ +1 /2 ( yH-V2|| 6 ||^|| B ||^ 



Summing over N yields 



m A" 1 



E E M NA )\ Y N < CE E Y^MNA^iNA^ + CA-^S-^WyUoo (3.15) 

i=l N=l i^j N=l 

+ CA-^S^WaUoo + CA^-^-^WbWsWBWg . (3.16) 

Note now that by Lemma l3~71 one can hope that Y^ m T for every i and N. Therefore, the left 
hand side is very close to T/A times the sum of the L 1 norms of the 6j. More precisely, one has 
for every i, the bound 



\aJ2Mna)\ - \\bi\\ L i 



N=l 



< INIhA 11 , 



so that, multiplying (13. 1 5 1 > by A/T, we get the bound 



E NUi < CA 6 \\b\\ s + CS- H \\y\\ L <~ + CA 6 6 6 - H \\b\\ 6 \\B\\ 6 
i=i 

m A' 1 

+ CS^Wah- + CA 1 ' 2 ^ 2 - 11 J2 E M NA Wn " T\ 



i=l N=l 



+ CAV^V2-h £ J- Y^\b t (NA)b 3 (NA)\ 

i^j N=l 

< CA M \\b\\ij + C6- H \\y\\ L °° + CA il 5 A ~ H \\b\\j i \\B\\j i 

m A" 1 

+ CS X ~ H \\a\\ L ~ + CAW/2-H|| &b J- E ft? " > 

ij=l AT=1 

where we used Sij to denote the Kronecker delta. At this point, we note that, for every 7 < 1, one 
has the interpolation inequality 

||&lk»<C( 7 ||6||H + 7~ 1/ ^ll%i) • (3-17) 
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Therefore 



+ C 7 - 1 ^AV-^||6||^||i?||^ + C7||6||^ 



m A 

+ C 7 - 1 /^ A i/2 (5 i/2^|| & ||^ £ IY^'-T^i 

ij=l jv=1 



We now make the following choices for 7, 8, and A: 

H(l-H) 1 1-g 

7 e(i+H)(2-H) , ,5 £ h(2-h) f A ps e^p-s) . 

Note that, provided e is small, one has indeed i5 « A « 1. With these choices, and by choosing 
H sufficiently close to H, we see that there exists a constant a > and a constant C > such 
that 

IHI* 1 + 11*11* + £ gV2 Ag - 3 /2 I + e_1 llyll^ + 11°^ 

< c E - (1 + M| + |*||} + f £ ) + + ||„| U . 

Actually, the constant a can be brought arbitrarily close to (^ + h){2-h) — — Note tnat 
Lemma 13721 yields the bound 

\Y^-T5 t3 \>h)<C^(-c J ^ wzi y 
for every possible value of i, j, and N. This immediately implies 

Therefore, there exists a constant c such that, for e small enough, one has the bound 

\b\\ L ~>e^ and ||y|k- < e ) < p( £ S "1^ > ^ ) ^ 19) 

+ F( \\b\\l + \\B\\l + \\a\\ L «> > ce~ a / 2 ) < C p e? . 



The last inequality is obtained by combining d3- 1 8I > with the a priori bounds on the processes a, 
b, and B. We now turn to the proof of (I3.13I) . Fix again some small value of A to be determined 
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later. We then have for every t 6 [0, 1 — A] the inequality 

r-t+A 

a(s) ds 



< 2\\y\\ Lao + \(b(t), B(t + A) - B(t))\ + C\\b\\ 6 \\B\\ 6 A 2H 

< 2\\y\\ Lao + ||%oo||S|| # A 5 + C||%||5||£A 2 " . 



It is easy to show that, similarly to d3.17t one has the inequality 

r-t+A 

a(s) ds 



|«||l°° < 2A 1 sup 

te[o,i-A] 



+ 2A^||a 



ii 



so that 



|a|| L » < 2A^ 1 || 2/ || L o + CWbW^WBWfjA^- 1 + Cllbll^llSl^A 2 ^- 1 + 2A & \\a\\ 6 



< C(A- 1 (||y|| ioc + 116111.) + A^OIBUI + + ||a||^)) . 



At this point, we choose A ps e a / H with a as in (13.191 . This implies (by choosing as before H 
sufficiently close to H) that there exists (3 > such that 



|a||i°° < C<r 



+ 

e e 



h\\ 2 

0\\ L oc 



+ 11*111+ 



Therefore, for e small enough, there exists a constant c such that 

H|l°° > e^ 2 and < e) < P( > e°^ 2 and ||y||z,°° < e 



+ 



1^111 + ' 



ifl 



for arbitrary values of p. The last inequality is obtained by combining d3.19t with the a priori 
bounds on the processes a, b, and B. This concludes the proof of d3.13t and thus of Proposition 
l3~4l □ 



4 Existence and smoothness of the density under Hormander's con- 
dition 

We now arrive to the heart of our study and are interested in the study of the existence and regu- 
larity for the density of solutions of stochastic differential equations on IR n 

Xf = x+ f* V (X*)ds + Vf Vi(X*)dBi (4.20) 
Jo i=1 Jo 

where the V^'s are C°° -bounded vector fields on M. n and B is the d dimensional fractional Brow- 
nian motion defined by (I2.lt . For this type of equations, existence and uniqueness of the solution 
have been investigated by many authors (for instance Nualart-Rascanu in [NR02 1). It has also been 
shown very recently in [ HN06 1 that the law of the solution to d4.20t possesses a smooth density 
with respect to the Lebesgue measure in the elliptic case. Recall first the following a priori bound 
on the solutions to d4.20b . (See for example BNR02I for a proof.) 
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Lemma 4.1 For every 7 < H and every p > 0, E(||X||^) < +00. 

Let us now denote by & the stochastic flow associated with equation d4.20t . that is <3?t(x) = Xf. 
From IHN06I . we can deduce: 

Lemma 4.2 The map <&t is C 1 an d the first variation process defined by 

T d$ t 
-fo^t — -5—, 
ox 

satisfies the following equation: 

Jo^t = W R n + f DV Q (Xf)3 ^ s ds + Vf DVi(Xf)3 ^ s dBi. 
Jo i=1 Jo 

and, for every p > 1, 

e(h ^ p) < +00. 

Furthermore, for every i = 1, . . . , n, t > 0, a?i<i x E M n , X^' 4 G B°°('H). Moreover, 

D>,X? = J ^tJol s Vj(X s ), j = l,...,d, 0<s<t, (4.21) 

where Y)iX^' 1 is the j-th component ofT) s X^' 1 . 

We can now turn to our version of Hormander's theorem for stochastic differential equa- 
tions that are driven by a fractional Brownian motion with Hurst parameter H > 1/2. If I = 
(ii, . . . , ik) £ {0, . . . , d} k , we denote by Vj the Lie commutator defined by 

V I = [V h ,[V i2 ,...,[V ik _ 1 ,V ik }...}. 

We also define the sets of vector fields 

n 

V n = {V T , Ie {0,...,d} x {l,...^}™" 1 } , v n = (J V fc . 

k=0 

With these notations, the main result of this article is the following: 

Theorem 4.3 Assume that, at some xo £ M n , there exists N such that 

span{l/(x ), V G V N } = R n . (4.22) 

Then, for any t > 0, the law of the random variable X^° has a smooth density with respect to the 
Lebesgue measure on W 1 . 

Given the results from the previous sections, the proof of this theorem is by now quite standard 
and follows closely the argument given for instance in BBisSll . BNor86l or INua95l . The main 
difference is that it is not a priori obvious how to relate the L°° bounds obtained in Proposition 13 .41 
to the fractional Sobolev norms appearing in the statement of Theorem 12. II 
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Lemma 4.4 Let if > 1/2 and let Ti be defined as above. Then, for every 7 > H — 1/2 there 
exists a constant C such that 

1 3+1/7 



H > L 



2+I/7 ' 

7 



/or every continuous function f £ Ti. Here, ||/|| 7 denotes as before the ^-Holder norm of f. 
Proof. Let denote the fractional derivative of order a of /, defined by 



1 d 







We also introduce the operator V°: defined by 



Vtf{t) = ~Y^T J (S - tr a ^(f(s) - f(t)) d8 , 

which is nothing but the adjoint of V a in L 2 (R + ). 

Since T a and V a are each other's inverse [SKM93|, ( 12.31 ) implies by Cauchy-Schwartz that 

\(f,g)v>\ = KX^- 1 / 2 /^-" 1 ^)^! < \\f\\n\\V H ^ l2 g\\ • (4-24) 



The problem with (14.241) is that we would like to apply it to a function / which is 7-Holder 

i<2 

continuous on [0,1], but does not necessarily vanish at either or 1, so that P_ / does in 
general not belong to L 2 . If we define however h(t) = t 7 (l — t) 7 for t G [0, 1] and h(t) = 
for t > 1, then ||//i|| 7 < C||/|| 7 , but fh vanishes at and at 1. In particular, this implies that 

\\V H ~ 1/2 fh\\ <C||/|| 7 ,sothat 

min - c — re — ' 

for some constant C. On the other hand, it is a straightforward calculation to check that 

1 3+1/7 

t y(\-tyf(t)dt> c 



1+1/7 ' 

7 



which implies the desired result. □ 

Corollary 4.5 Let y be a random process with sample paths that are almost surely ^-Holder 
continuous for some 7 > H — 1/2. Then, there exists an exponent a such that 

W{\\y\\n <e)< P(||y||z~ < e a ) +P(|MI 7 > e~ a ) , 
for every e sufficiently small. 
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Proof. It follows from Lemma l4~4l that 

P(||y|| w < e )<F(C|| y ||^|| y ||^ 1 /7 <e) . 

For an arbitrary pair of positive random variables X and Y, one always has 

F(X/Y <e)< F(X < e l ~ a ) + F(Y > e~ a ) , 

so that the claim follows by taking a small enough. □ 

This provides us with the necessary tools to complete the 

Proof of Theorem \4.3\ We shall show that X x ° admits a smooth density with respect to the Lebesgue 
measure, by using the Malliavin covariance matrix Fi associated with Xf°. Note that we can con- 
sider the case t = 1 without any loss of generality by rescaling the vector fields Vi appropriately. 
Let Ti be the Malliavin covariance matrix associated with X x °. By definition, we have 



ri = ((Dxi' xo ,r>x{' xo ) n ) 



l<i ,j<n' 



To show that X®° has a a smooth density, it suffices to show that with probability one Ti is 
invertible and that for every p > 1, 

E( — - , ) < +00. 
VI detTi \p J 

From Lemma l4~2l 

BiXf = Jo^J^pQ, j = l,...,d, < 8 < 1. 

Therefore, 

r x = H(2H - 1) Jo^ T !\ l u V{X^)V(X^) T {^l v ) T \u-v\ 2H - 2 dudvJ^, 
Jo J 

where V denotes the n x d matrix {V\ . . . V4). Since Jo^i is almost surely invertible with inverse 
in LP, p > 1, in order to show that Ti is invertible with probability one, it is enough to check that 
with probability one, the matrix 

Ci = f t J iuV(K )V(X x v o) T (3-i v ) T \ u-v \ 2H ~ 2 dudv 
Jo Jo 

is invertible and satisfies for every p > 1, 

Ki^kf) < + °°- (4 - 25) 
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By using Proposition 13.41 the idea is now to control the smallest eigenvalue of C\ by showing 
that it can not be too small. More precisely, recall (Lemma 2.3.1. in BNua95IO that if for any p > 2, 
there exists eo(p) such that for every e < eo(p), 

sup F((v,C lV ) <e)< e p , 

then C\ is invertible with probability one and (I4.25t holds for every p > 1. We thus want to 
estimate F((v, C\v) < e). Let us observe that 




(v, (& a Vi){x ))(v, mVjXxo)) dsdt 



= £ll<«.(^)(*o))l&. 
J'=l 

Fix now an (arbitrarily large) value p > 0. It follows from Lemma PTTl and Corollary I4.5l that there 
exists a > such that 

H(v,C lV } <e)< Ce p + min P(\\(v,(^*Vi)(x ))\\ LOO < e a ) . (4.26) 

i=l,...,d 

Note now that if V is an arbitrary bounded vector field with bounded derivatives, the chain rule 
reads 

mv)(*o)= [\y,mvo,v))(xo))ds+j2 f{ y ,mvvV])(^))dBi. 

Jo j=1 Jo 

It thus follows from Proposition 13 .41 and Corollary I4.5l that there exists a such that 

F(\\(v,(<S>*V)(x ))\\ L oo < e ) <CeV+ min P(||( W) ($* [V h V])(x ))\\ L °o < e a ) (4.27) 

i=0,...,d 

Consider now the integer N from the assumption. Combining (14.26t with ( 14.271 . we see that there 
exists a > such that 

H{v,C lV ) <e) <Ce p + min F(\\(v, ($* V)(x ))\\ Loo < e a ) . 

v&v N 

for all e small enough. On the other hand, we know by assumption that {V{xq) , V G Vn} 
spans all of W 1 , so that there exists some V G Vat such that («, V(xo)) 7^ 0. Therefore, one has 
F((v , Cif ) < e) < Ce p for all e sufficiently small, which is the required bound. □ 

5 Asymptotics of the density in small times 

In order to obtain asymptotics of the density of hypoelliptic diffusions on the diagonal in small 
times, one method consists to approximate the diffusion by the lift of the Brownian motion in a 
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nilpotent Lie group that is called a Carnot group (see IBau0 5 ] or IB A89I ). In a recent work IBC06I . 
Baudoin and Coutin have introduced and studied fractional Brownian motions on Carnot groups. 
We shall see that the solution of a stochastic differential equation driven by fractional Brownian 
motions can, under Hormander's type assumptions, be approximated by fractional Brownian mo- 
tions on Carnot groups. From this approximation, we will deduce an asymptotic development of 
the density in small times. 

We recall first the notion of Carnot group (see e.g. IBau05IO and the main results that are 
obtained in IBC06I . 

Definition 5.1 A Carnot group of step (or depth) N is a simply connected Lie group G whose Lie 
algebra can be decomposed as 

Vi ffi . . . V N , 

where 

[V i ,V j ]=V i+j 

and 

V s = 0, for s > N. 

Notice that the vector space V\, which is called the basis of G, Lie generates g, where g denotes 
the Lie algebra of G. It is possible to show that for every N > 1, up to an isomorphism, there 
exists exactly one A^-step nilpotent Carnot group with basis M. d . This group shall be denoted by 

G N (R d ). 

Since a Carnot group G is A^-step nilpotent and simply connected, the exponential map is 
a diffeomorphism. On g we can consider the family of linear operators St : — > 0, t > 
which act by scalar multiplication t % on Vi. These operators are Lie algebra automorphisms due 
to the grading. The maps St induce Lie group automorphisms At : G —> G which are called the 
canonical dilations of G. Let us now take a basis U\, ... ,Ud of the vector space Vi. The vectors 
Ui can be seen as left invariant vector fields on G so that we can consider the following stochastic 
differential equation on G: 

d ,-t 

dX t = Y. Ui{X s )dBl t > 0, (5.28) 
i=i J ° 

which is easily seen to have a unique solution associated with the initial condition Xq = 1q. The 
driving process (B t )t>o is here a fractional Brownian motion with Hurst parameter H > 1/2. The 
process (X t ) t >o is called the lift of (B t )t>o in the group G. For this equation, the assumptions of 
Theorem l4.3l are obviously satisfied, so that we have a smooth density for X t ,t > 0, with respect 
to the Haar measure of G. 

We have then the global scaling property 

(X c t)t>o = (A c ifXt)t>o- 

(See IBC06I for a proof.) This scaling property leads directly to the following value at 1^ of the 
density p t of X t with respect to the Haar measure of G: 

C 

PtO-G) = TW * > °> 
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where C > and D = Yli=i * dim From this, we will see how to deduce asymptotics in small 
times for any hypoelliptic stochastic differential equation driven by fractional Brownian motions. 

From now on, we consider d vector fields Vi : R n — > R n which are C°° bounded and shall 
always assume that the following assumption is satisfied. 
Strong Hormander's Condition: For every x S R n , we have: 

span{VK*U G U fc >i{l, . . . , d} k } = R n . 

We recall that if I = (ii, . . . , i^) G {1, . . . , d} k is a word, we denote by Vj the commutator 
defined by 

V I = [V h ,[V i2 ,...,[V ik _ 1 ,V ik ]...]. 

Let us introduce some concepts of differential geometry. The set of linear combinations with 
smooth coefficients of the vector fields Vy , . . . , is called the differential system (or sheaf) gen- 
erated by these vector fields. It shall be denoted by V in the sequel. Notice that V is naturally 
endowed with a Coo(R n , R) -module structure. For x S R n , we denote 

v(x) = {x(x),x e v}. 

If the integer dimP(x) does not depend on x, then V is said to be a distribution. The Lie brackets 
of vector fields in V generate a flag of differential systems, 

V = v l CD 2 c ■■■ <zv k <Z ■■■ , 

where V k is recursively defined by the formula 

v k = V k ~ x + [V,V k ~ 1 }. 

As a module, V k is generated by the set of vector fields Vi, where / describes the set of words 
with length k. Moreover, due to Jacobi identity, we have [V 1 , P J ] C V l+J '. This flag is called 
the canonical flag associated with the differential system V. Hormander's strong condition, which 
we supposed to hold, states that for each x € R n , there is a smallest integer r(x) such that 
-pr(x) _ p Qr eac fo x g jg,n^ t ^ e canon j ca i fl a g induces a flag of vector subspaces, 

V{x) C V 2 (x) C • • • C V r(x \x) = R n . 

The integer list (dim P fc (x)) 1<k<r ^ is called the growth vector of V at x. The point x is said to 
be a regular point of V if the growth vector is constant in a neighbourhood of x. Otherwise, we 
say that x is a singular point. On a Carnot group, due to the homogeneity, all points are regular. 
Let Vi = V 1 /V 1 " 1 denote the quotient differential systems, and define 

M{v) = Vi © • • • e v k © • • • . 

The Lie bracket of vector fields induces a bilinear map on M{V) which respects the grading: 
[Vj, Vj] C Vj+j. Actually, N(T>) inherits the structure of a sheaf of Lie algebras. Moreover, if x 
is a regular point of V, then this bracket induces a r(j;)-step nilpotent graded Lie algebra structure 
on J\f(V)(x). Observe that the dimension of J\f(V)(x) is equal to n and that from the definition, 

(Vi(x), ...,V d (x)) Lie generates M{V)(x). 
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Definition 5.2 If x is a regular point of V, the r(x)-step nilpotent graded Lie algebra N(V)(x) 
is called the nilpotentisation of V at x. This Lie algebra is the Lie algebra of a unique Carnot 
group which shall be denoted Gr(V)(x) and called the tangent space to V at x. The integer 
D = X)I=i kdimV k (x) is called the homogeneous dimension of Gr(V)(x). 

Proposition 5.3 Let x be a regular point of D. Let pt, t > 0, denote the density with respect to 
the Lebesgue measure of the solution of the stochastic differential equation 

X? = x + Y,f Vi{X x s )dBl (5.29) 



where (Bt)t>o is a d-dimensional fractional Brownian motion with Hurst parameter H > 1/2. 
We have, 

where C(x) is a non negative constant and D{x) the homogeneous dimension of the tangent space 
Gr(V)(x). 

Proof. Let us first introduce some notations: For k > 1, we denote by A fc [0,t] the simplex of 
ordered /c-tuples with values in [0, t], i.e. 

A k [0,t] = {(t!,...,t k ) G [0,t] k ,h < ... <t k }; 

If I = . . . ,ifc) G {1, . . . , d} k is a word with length k, we define the corresponding iterated 
integral of B by 

f dB ! =[ dB£---dB%, 

JA k [0,t] J0<ti<...<t k <t 

where the right hand side consists of nested Riemann-Stieltjes integrals. 

We denote &t the group of the permutations of the index set {1, . . . , k} and if a G and I 
is a word I = (ii, . . . , i k ), we denote by a ■ I the word (i a ru , . . . , i a (k))- 

If a £ 6fc, we denote e(a) the cardinality of the set 

{j e{l,...,k-l},a(j) >a(j + l)}, 

i.e. e(a) is the number of raising sequences of a. Finally, if I = (ix, . . . , i k ) G {1, . . . , d} k is a 
word 

Ne(cr) 



e(a) 



dB a 11 . 



As a consequence of Proposition 23 in IFV06L we get the following approximation result 



X? 



r(x) 

ex p ( E E 

k=l I=(i x ,...,i k ) 



(x) +t~~5~R(t), i > 
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Here, the remainder term R(i) is such that there exist a, c > such that, for all A > c, 
limsuppf sup Z^ 1 ' | Rat( S ) |> A^ +1 ) j < exp — V 



It is shown in IBau051 Proposition 3.7, pp 66] that one can write 

{x)=F(Bt) t 



(r(x) 
2 £ Aj(S) t Vj 
fe=l I=(ii,...,i k ) 

where {B^) t >Q is the lift of (B t )t>o in the free Carnot group G r (,,.)(R d ) and where F is a map 

G r ( a: )(lR ,i ) —* K n 

that is Pansu differentiable at Ik , , . Furthermore, its Pansu's derivative 

d P F(l Gr(;c) ) : G r{x) (R d ) - Gr(2?)(x) 
is a surjective Carnot group morphism. Now, we have 



and 



^ « xA^ V ^d P F(0 Gr{x)(Rd) )(Bl), 
in small times. Since dpF(0 Gr is surjective, the random variable 

admits a density with respect to the Lebesgue measure on M. n . This concludes the proof. □ 
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